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I.  Introduction 


Sound  propagation  in  the  atmosphere  is  affected  by  vertical  temperature 
and  wind  gradients.  These  gradients  produce  a  refracting  medium  in  which 
sound  speed  is  a  function  of  height.  A  single  sound  speed  profile  as  a 
function  of  height  above  the  ground  can  be  obtained  from  individual  temp¬ 
erature  and  wind  gradients  via  empirically  derived  equations  or  measured 
directly  or  indirectly  with  sounders.  The  temperature  and  wind  effects  can, 
therefore,  be  considered  collectively,  for  acoustic  purposes,  as  a  variation 
of  sound  speed  with  height. 

A  mathematical  sound  ranging  model  must  effectively  Incorporate  vertical 
sound  velocity  gradients.  The  computer  programs  developed  to  model  sound 
propagation  using  ray  techniques  presented  by  Roth'*'  require  as  input  either 
the  vertical  temperature  and  wind  profiles  or  the  vertical  sound  velocity 
profile. 

It  is  assumed  that  current  meteorological  vertical  temperature  and  wind 

velocity  profiles  will  be  utilized  in  the  sound  propagation  model.  However, 
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when  current  meteorological  data  are  unavailable,  models  that  have  been 
developed  to  extend  limited  surface  temperature  and  wind  conditions  vertically 
may  be  of  value.  This  report  presents  several  temperature  and  wind  profile 
models  and  suggestions  for  their  use. 

II.  Wind  Models 
1.  Power  Wind  Law 

The  power  law  is  often  used  by  engineers  to  portray  the  vertical  wind 
profile  (p-profile  theory) : 


where  is  the  wind  speed  at  the  anemometer  height  Z^,  and  Uj  is  the  predicted 
wind  speed  at  height  Z The  power  p  is  taken  to  be  constant  with  height  in 


a  given  range.  Different  values  of  p  should  be  chosen  depending  on  the  height 
range  over  which  the  power  law  is  applied. 


A  useful  method  for  estimating  p  can  be  derived  from: 

p  -  1/ln  (Z/Zo) 

where,  in  practice,  one  chooses  Z  as  the  geometric  mean  height  in  the  layer 

to  which  the  power  law  profile  is  to  be  applied.  Zo  is  the  roughness  length 
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which,  according  to  Wagner  varies  from  1x10  m  for  a  mud  flat  to  2.83m  for 
a  fir  forest.  For  example,  in  smooth  terrain  (Zo«0.01m)  in  a  layer  of  mean 
height  xOm,  p  is  about  1/7,  a  value  commonly  suggested  in  engineering  books. 

But,  over  rough  terrain  (Zo“lm,  Za10m) ,  p  is  much  larger. 

a 

According  to  Hans  A.  Panofsky  ,  the  p  -  profile  theory  for  wind  profiles 
does  not  fit  very  well  when  the  terrain  Is  uniform  and  mechanical  turbulence 
prevales.  Even  then,  a  power  law  can  be  fitted  reasonably  well  over  a  small 
height  range.  The  use  of  p  -  theory  does  not  permit  variation  of  wind 
direction  with  height. 

2.  Logarithmic  wind  law 

Like  the  power  law,  the  logarithmic  wind  law  for  modeling  the  vertical 
structure  of  the  wind  applies  only  in  the  lowest  turbulent  layer.  Hess'* 
defines  the  logarithmic  wind  law: 

TT  U*  ,  Z 

U  “  —  In  — 
k  Zo 

where  k  is  a  pure  number  called  the  Von  Karmen  constant  and  Zo  is  the  roughness 
parameter.  Von  Karmen's  constant  is  nondlmensional  and  has  been  found  to  be 
approximately  0.4.  U*,  the  friction  velocity,  may  be  approximated,  according 


to  Blackadar  ,  by: 


U*  •  k  U 


SFC 
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where  Ugp^  is  the  surface  wind  speed  and  Za  is  the  height  of  the  anemometer. 
When  the  lapse  rate  is  not  neutral,  the  logarithmic  wind  law  is  not  a 
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good  description  of  the  wind  profile.  However,  even  through  the  real  atmos¬ 
phere  is  seldom  neutral,  it  has  been  demonstrated  that  the  logarithmic  wind 
law  tends  to  give  more  realistic  results  than  the  power  law  in  Che  lowest 
turbulent  layer  which  is  normally  less  than  100m  thick. 

3.  Ekman  Layer  Model 

Above  50  to  100m  the  assumption  that  the  turbulent  stress  is  independent 
of  height  is  no  longer  valid.  Observations  demonstrate  that  the  stress  de¬ 
creases  with  height  until  it  becomes  negligible  and  the  wind  responses  primarily 
to  the  pressure-gradient  and  the  corlolis  force.  This  height  is  called  the 
gradient  level  and  is  found  near  1000m. 

Hess'’  derives  the  equations : 

-aZ 

u  ■  Ig(l-e  cos  aZ) 

-aZ 

v  ■  Ug  e  sin  aZ 

where  Ug  is  the  geostrophic  wind  speed  and  a  •'/  f/2K. 
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K  is  a  constant  and  is  approximately  5  x  10  cm  sec  in  the  middle 
latitudes. 

The  consequence  of  these  equations  is  that  they  demonstrate  that  as  the 
wind  turns  in  response  to  the  trigonometric  terms  of  the  equations,  the  amplitude 
of  these  terms  decrease  exponentially.  Thus  the  wind  speed  and  direction 
converge  rapidly  on  the  geostrophic  values . 

III.  Vertical  Temperature  and  Wind  Profile  Model 


A  number  of  models  exist  to  compute  the  vertical  distribution  of  temperature 
and  wind  with  height.  Reference  2  lists  a  number  of  vertical  wind  profile 
models  available  at  ASL.  Blackadar^  presents  another  model  which  incorporates 
a  routine  to  compute  the  vertical  distribution  of  wind  with  height  suitable  to 
a  neutral  temperature  profile  with  the  same  potential  temperature  as  the  surface 
temperature.  The  Blackadar  routine  was  modified  to  include  the  power  wind 
law  as  an  option.  Inputs  to  the  model  include  the  temperature  lapse  rate 
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(°C/m),  average  height  of  obstructions  (m) ,  surface  wind  speed  (m/s)  at  some 
height  (m) ,  latitude  (DEG),  and  the  vector  difference  between  the  surface 
wind  and  the  thermal  wind  (m/s). 

The  model  approximates  the  roughness  length  by: 

Zo  %  AVEHGT/30.0 

where  AVEHGT  is  the  average  height  of  the  obstructions  in  meters.  The  friction 
velocity  is  approximated  by: 

U*  -  0.4  *  SFCWS  [In  (ANEM0/Zo) ] -1 

where  SFCWS  is  the  surface  wind  speed  at  the  anemometer  height ,  (ANEMO) . 

For  barothropic  conditions,  the  vector  difference  between  the  surface 
wind  and  the  thermal  wind  is  zero.  The  X  component  of  the  vector  is  taken 
to  be  in  the  same  direction  as  the  surface  wind.  The  angle  between  the  surface 
wind  vector  and  the  thermal  wind  vector  will  be  defined  as  B.  A  surface  meteor¬ 
ological  chart  or  forecast  of  thermal  winds  for  the  time  and  place  in  question 
is  assumed  to  be  available.  In  many  instances  the  atmosphere  is 

essentially  barothropic  and  in  those  cases  only  small  errors  will  be  intro¬ 
duced  if  VDIF,  the  vector  difference,  is  set  to  zero. 

The  Blackadar  model  utilizes  the  logarithmic  wind  law  in  the  friction 
layer  and  the  Ekman  layer  model  in  the  planetary  boundary  layer.  The  pro¬ 
posed  model  predicts  temperature  and  wind  at  50m  intervals  up  to  1000m  which 
is  the  approximate  height  of  the  planetary  boundary  layer  in  mid- latitudes. 

IV.  Proposed  Use  of  the  Temperature  and  Wind  Profile  Model 

The  temperature  at  the  surface  and  temperature  lapse  rate  are  inputs  to 
the  model  for  determining  the  vertical  temperature  distribution.  The  default 
condition  on  the  temperature  lapse  rate  is  0.0065°C/m. 

The  surface  wind,  roughness  length,  friction  velocity,  thermal  winds,  and 
latitude  are  input  to  the  model  for  determining  the  vertical  distribution  of 
wind  speed,  u  and  v  component,  with  height. 
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For  use  of  the  power  law  model,  only  the  temperature  at  the  surface, 
temperature  laps  rate,  surface  wind  speed,  height  of  the  anemometer,  and 
average  height  of  the  obstructions  need  be  entered.  The  selection  switch 
(NSWTCH)  variable  must  be  non-zero  if  this  option  is  selected. 

Table  1  and  Table  2  provide  examples  of  the  subroutine  output  given 
realistic  values  for  surface  wind  (4.76  mile/hr  and  9.51  mile/hr)  and 
temperature  (10°C).  Given  the  specified  initial  parameters,  which  included 
VDIF  and  B,  a  driver  program  called  the  proposed  subroutine  to  compute 
the  vertical  temperature  and  vertical  wind  profile  using  the  power  law 
and  the  modified  Blackadar  routine. 

The  power  law  and  modified  Blackadar  routines  calculated  reasonable 
values  for  the  variation  of  wind  with  height.  The  Blackadar  routine  did 
produce  values  for  the  wind  that  were  about  twice  as  large  as  the  power 
law  results  at  1000  m  for  the  cases  tested.  Since  Blackadar  utilized 
equations  which  have  trignometric  terms,  the  resultant  wind  does  turn 
with  height. 

The  inclusion  of  the  vector  difference  term  demonstrates  the  effect 
it  has  on  the  model  output.  The  test  case  of  30°  shows  a  slight  increase 
in  the  U  component  of  the  wind  at  all  levels  compared  to  the  corresponding 
barotrophic  level. 

Appendix  1  provides  a  computer  listing  of  the  proposed  model  and 
Appendix  2  defines  the  model’s  input  variables. 
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TABLE  1 


Subroutine  Test  Results  Given  a  Friction 
Velocity  of  0.25  and  Roughness  Parameter  of  0.1 


Input  Parameters: 

OLAT  -  40.0 
TEMP  -  10.0 
DELT  -  0.0065 


SFCWS  -  2.126  M/S 
ANEMO  -  3.0  M 
AVEHGT  -  3.0  M 


WINDS 

(M/S) 

POWER 

LAW 

MODIFIED 

BLACKADAR 

VDIF  - 

0.0 

VDIF  - 

1.5 

B  - 

0.0 

B  - 

30.0 

HEIGHT 

TEMP 

U 

V 

u 

V 

U 

V 

(M> 

(°C) 

(M/S) 

(M/S) 

(M/S) 

(M/S) 

(M/S) 

(M/S) 

25.0 

9.84 

2.67633 

0.00000 

3.45132 

0.00000 

3.45132 

0.00000 

75.0 

9.51 

3.01538 

0.00000 

4.48661 

-0.25090 

4.48827 

-0.24658 

125.0 

9.19 

3.18734 

0.00000 

5.29716 

-0.34377 

5.30482 

-0.32687 

175.0 

8.86 

3.30594 

0.00000 

5.90503 

-0.50770 

5.92371 

-0.47193 

225.0 

8.54 

3.39739 

0.00000 

6.33695 

-0.71951 

6.37171 

-0.66035 

275.0 

8.21 

3.47222 

0.00000 

6.62257 

-0.95843 

6.67819 

-0.87288 

325.0 

7.89 

3.53577 

0.00000 

6.79149 

-1.20702 

6.87224 

-1.09334 

375.0 

7.56 

3.59114 

0.00000 

6.78131 

-1.45162 

6.98081 

-1.30900 

425.0 

7.24 

3.64027 

0.00000 

6.88640 

-1.68223 

7.02763 

-1.51058 

475.0 

6.91 

3.68449 

0.00000 

6.85736 

-1.89219 

7.03264 

-1.69191 

525.0 

6.59 

3.72475 

0.00000 

6.80089 

-2.07766 

7.01196 

-1.84948 

575.0 

6.26 

3.76173 

0.00000 

6.72995 

-2.23714 

6.97805 

-1.98194 

625.0 

5.94 

3.79593 

0.00000 

6.65415 

-2.37083 

6.94010 

-2.08955 

675.0 

5.61 

3.82779 

0.00000 

6.58024 

-2.48021 

6.90450 

-2.17377 

725.0 

5.29 

3.85760 

0.00000 

6.51255 

-2.56754 

6.87535 

-2.23680 

775.0 

4.96 

3.88564 

0.00000 

6.45358 

-2.63554 

6.85496 

-2.28125 

825.0 

4.64 

3.91210 

0.00000 

6.40438 

-2.68710 

6.84426 

-2.30991 

875.0 

4.31 

3.93717 

0.00000 

6.36502 

-2.72507 

6.84320 

-2.32549 

925.0 

3.99 

3.96100 

0.00000 

6.33485 

-2.75209 

6.85115 

-2.33054 

975.0 

3.66 

3.98370 

0.00000 

6.31284 

-2.77056 

6.86704 

-2.32735 
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TABLE  2 


Subroutine  Test  Results  Given  a  Friction 
Velocity  of  0.5  and  Roughness  Parameter  of  0.1 


Input  Parameters: 

OLAT  -  40.0  DEG 

TEMP  =  10.0  °C 

DELT  -  0.0065  °C/M 


SFCWS  =  4.251  M/S 
ANRMU  -  3.0  It 
AVEHGT  «  3.0  M 


WIND 

POWER 

LAW 

MODIFIED 

BLACKADAR 

VDIF  - 

0.0 

VDIF  - 

1.5 

B  * 

0.0 

B  = 

30.0 

HEIGHT 

TEMP 

U 

V 

u 

V 

U 

V 

(M) 

(°C) 

(M/S) 

(M/S) 

(M/S) 

(M/S) 

(M/S) 

(M/S) 

25.0 

9.84 

5.35140 

0.00000 

6.90101 

0.00000 

6.90101 

0.00000 

75.0 

9.51 

6.02935 

0.00000 

8. 33166 

-0.48419 

8.33179 

-0.48378 

125.0 

9.19 

6.37319 

0.00000 

9.36418 

-0.48306 

9.36592 

-0.47828 

175.0 

8.86 

6.61032 

0.00000 

10.28351 

-0.53223 

10.28893 

-0.51869 

225.0 

8.54 

6.79318 

0.00000 

11.09333 

-0.62590 

11.10474 

-0.59970 

275.0 

8.21 

6.94280 

0.00000 

11.79865 

-0.75807 

11.81853 

-0.71585 

325.0 

7.89 

7.06988 

0.00000 

12.40554 

-0.92277 

12.43644 

-0.86162 

375.0 

7.56 

7.18058 

0.00000 

12.92078 

-1.11415 

12.96525 

-1.03161 

425.0 

7.24 

7.27883 

0.00000 

13.35161 

-1.32665 

13.41217 

-1.22067 

475.0 

6.91 

7.36726 

0.00000 

13.70552 

-1.55505 

13.78462 

-1.42396 

525.0 

6.59 

7.44776 

0.00000 

13.99005 

-1.79453 

14.09004 

-1.63702 

575.0 

6.26 

7.52169 

0.00000 

14.21266 

-2.04072 

14.33573 

-1.85578 

625.0 

5.94 

7.59008 

0.00000 

14.38057 

-2.28975 

14.52880 

-2.07663 

675.0 

5.61 

7.65377 

0.00000 

14.50072 

-2.53818 

14.67600 

-2.29639 

725.0 

5.29 

7.71338 

0.00000 

14.57965 

-2.78308 

14.78371 

-2.51233 

775.0 

4.96 

7.76944 

0.00000 

14.62346 

-3.02195 

14.85788 

-2.72212 

825.0 

4.64 

7.82236 

0.00000 

14.63781 

-3.25273 

14.90397 

-2.92386 

875.0 

4.31 

7.87250 

0.00000 

14.62783 

-3.47378 

14.92700 

-3.11602 

925.0 

3.99 

7.92014 

0.00000 

14.59819 

-3.68380 

14.93143 

-3.29741 

975.0 

3.66 

7.96553 

0.00000 

14.55307 

-3.88186 

14.92132 

-3.46717 
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* 

5LeR0UTlN£  wCTEMP  USl»TCH,CLAT,TEMP,CEcT,SFCwS,.AN&M0,.AvEhG7»VDIF,e) 

PROGRAM  tc  calculate  temperature  anc  »INC  profiles 

USING  A  MCCIPIEC  VERSKN'CF  BLACKCAP'S  ROUTINE 


EKAN  LAYER  IMCCEL  TC  CALCULATE  taINC  FRCFILt  GIVEN  FRICTICN 
VELOCITY,  IROCGKNESS  LENGTH,  LATITUC6,  SURFACE  TEMPERATURE • 
SURFACE  *»INC»  AND  .SCALED  HEIGHT  CP  THfc  SURFACE  LATER 

QLAT  (DEG)  LATIHLE  OF  AKEA 

AVEHGT  ( M )  A  VERAGE  *  HE IGHT  >CF  CbSTRUCTlQNS 

SF  CWb  (M/S)  SURFACE  tltrC  SPEED 

ANEMO  <M>  iHEIGHT  CF  THE  ANEMCMtlEM 

VC  IF  (M/S)  SURFACE  UNC  VECTDR  MINUS  THE  THERMAL 

B  (DEG)  ANGLE  BETWEEN  TiHE  SURFACE  Wlfi'D  VECTOR 

In  INC  VECTOR 

TEMP  t  C)  SURFACE  TEMPERATURE 

DELT  TEMPERATURE  lapse  rate 


»  UNO 

And  THE  THERMAL! 


GCMMUN  D  CEO)  »V  (2  0)  »T  UC)  »2HT<20>  *ZtC  (40)  ,2  (SO) 

BE  =  SCAL&C  iHElC-HT  .CF  THE  SURFACE  LA YUR 
BEsO.Ol 

BETA  s  D , V  0CMFCN6NT  'CF  T.HE  UH-UG  VECiCR 
BETAJtn.lfi 

VK*0 .4 
VK2=VK**2,C 
l  02*50. 

ZTH*U2/2.C 
ZH  T  ( I J  *ZTH 
B=e/57,2se 
SX=0,0 
SYsO.O 

Z2ERUsAVfaHGt/3C.O 
USTRs0.h*SF0«S/ALOG (ANEMC/ZZERC) 

IF (VUlF.tG,C,0)  GO  TC  2 
S/sVCIF«CGS  (-fc  ) 

SY*VUIF<*SI)N  {.fc  > 

\c  IF  (DtLT  ,NE  *0,0)  GO  TC  3 
DELTsO.0065 
13  CCNTINUE 

IF (USTR.LE.0,0)  GC  TC  195 

WRITt (6,11 )CU AT ,USTR,Z2ERC, TEMP, SA,«Y 

F=1,*14E-C4  '*  S  IN  (ICLAT/S7.29G) 

ZEC  (1)s(ZTih*F)/lST:R 
6C  DC  6b  1=2,20 

2HT  (i)*4MT  (I-D  *02 
ZECU)*UM  (I)«F)/USTR 
Z ( 1 ) 4UZ/2  *C 
Z(I)=z(I-l)yU2 
65  CCNTINUE 

IF (KSwTCH jNE»0)  GO  TC  5 
9C  OC  9b  1=1,20 

TEMPERATURE  .CALCULATES  FRCM  INPUT  CAFSE  RATE 

OEFAClT  CCNL 1 T 1CN  IS  THE  CRY  .ACIAOaTIC  LAHBt  RATE  O.C065  (C/M ) 
T(1)=TEMP-CEUT«2RT  (I) 

IF  (ZfeD(I)  *GT.0,01)  CC  TO  '51 

within  surface  layer  calculate  lcoritmmic  profile 
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U<I)  =  <USTR/VI‘  )«ALOC  (2ihT  { 1 )  /ZZEfiQ) 

V(I)=0,0 
SC  TO  95 
91  CONTINUE 

SX.ST  ARE  if/M  THE  THERMAL  «« iiNC  C*  Tht  6E0|T80BhIC  KlNC  -ShEAR 

OCE  TO  BABCCLlMCm 

U0  =  -o£TA^  (.5x«S>)/  CE.0«6ETA*VK2) 

VC*  otTA»  (.Sx-S>  }/  U.(l*b£TA*VK2) 

TERM!* (l.C/VN  )«AlOC  <C.01*LSTti/  (F*<UEPO)  ) 

1  -JL  C  ♦  ISA/VK2  )'•  (2tU  U )  -tt) 

TERMes  tAR  t-tETA*  t2EC  U>  -EE)  ) 

TEPHJs (UO«COS (EETA«  UtC II)-BE> ) 

1*V0*5IN<&ETA«  <2£D (!)-*£)>> 

U(I  )?(TthM*TEBh2«TtHK3)»LST8 
TEFMl  s  -VC* (SY/V^Z ) * <2EG (I)-BE) 

TEFMj*(VO-'COE  (E£TA«  UEC  ( I ) -BE  >  )  -10 
IAS  IN  (BETA*  (Zfcc  (l)-'fcE  )  >  > 

V  ( I )  9  <  TERM  !♦  1  ERM2*Tt  RF.3)  «LST8 
95  GCNT1NUE 
9fc  DC  1UO  1*1. 2C 

WPITt  (6.99)  1  ,2ihT  (I)  »T  (I)  »U ( I) .V  (1) 

IOC  CCNTlNOE 
RETURN 

.5  P=1.0/ALCe  (1C0G.0/22EK) 

DC  6  1*1.20 
T{I)*T£hP-C£LT*2hT (1) 

U(I)sSFC*S«  Khl  (I)/ANkhU)**P  . 

V<I)5O.0 
•6  OCNTiNCiE 
GC  TO  96 

11  PCRMaT  <  »0  *  . »LAT*  *  *F 10.4.5A. »tS JR*  *  *F 10 .A  »5X» *22ERCs  t»F10,A.5Xi 

1 • TEMP* *  »F 10  * A  »2 X  *  *  5*s » . F 1 0 .4 , 5* * *S Y» * »F 1 0 .*) 

99  FORMAT  ( •  •♦Ii:-»F10.1»FiC.2.2Fl0.5) 

199  STOP 

bnc 
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MODEL  INPUT  PARAMETERS 


SYMBOL 

DIMENSION 

DESCRIPTION 

OLAT 

DEG 

Latitude 

TEMP 

°C 

Surface  temperature 

DELT 

°C/m 

Temperature  lapse  rate 

SFCWs 

m/s 

Surface  wind  speed 

ANEMO 

m 

Height  of  anemometer 

AVEHGT 

m 

Average  height  of  obstructions 

VDIF 

m/s 

Vector  difference  of  wind  speed 
minus  thermal  wind 

B 

DEG 

Angle  between  surface  wind  and 
thermal  wind 

KSWTCH 

Model  selector  switch 
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t 


